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We give an account of the physical behaviour of a quasiparticle horizon due to non- 
Lorentz invariant modifications of the effective space-time experienced by the quasipar- 
ticles ("matter") for high momenta. By introducing a "relativistic" conserved energy- 
momentum tensor, we derive quasi-equilibrium states of the fluid across the "Landau" 
quasiparticle horizon at temperatures well above the quantum Hawking temperature. 
Nonlinear dispersion of the quasiparticle energy spectrum is instrumental for quasipar- 
ticle communication and exchange across the horizon. It is responsible for the establish- 
ment of the local thermal equilibrium across the horizon with the Tolman temperature 
being inhomogeneous behind the horizon. The inhomogeneity causes relaxation of the 
quasi-equilibrium states due to scattering of thermal quasiparticles, which finally leads 
to a shrinking black hole horizon. This process serves as the classical thermal counter- 
part of the quantum effect of Hawking radiation and will allow for an observation of 
the properties of the horizon at temperatures well above the Hawking temperature. We 
discuss the thermal entropy related to the horizon. We find that only the first nonlinear 
correction to the energy spectrum is important for the thermal properties of the horizon. 
They are fully determined by an energy of order -Epi an ck(T/-Epianck) 1/ ' 3 , which is well 
below the Planck energy scale -Epi anc k, so that Planck scale physics is not involved in 
determining thermal quantities related to the horizon. 

1. Introduction 

The quantum phenomenon of Hawking radiationEI is of kinematical origin, and 
relies on the existence of a Lorentzian signature background metric for the prop- 
agating particle, i.e., it depends upon the causal structure of space-time. It does 
not depend on the field equations determining the metricB and is, in particular, in- 
variant under conformal transformations of the metricB. Lorcntz invariance, which 
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is a symmetry of nature experienced to hold on large enough scales, may not be an 
exact symmetry on the smallest scales. Because a particle at an infinite distance 
from the black hole horizon experiences an infinite blueshift when traced back to 
the horizon, the horizon is a suitable means to probe the high frequency structure 
of particle spectra: Black hole physics is necessarily being influenced by the large 
momentum modes of the quantum fields. 

To trace the influence of the ultraviolet region upon the Hawking radiation from 
the black hole horizon, the condensed matter analogue of a black hole provides an 
orientation guide, with known underlying physical lawsBB. In condensed matter, 
ultraviolet deviations from "Lorentz invariancc" naturally arise, since the Lorentz 
symmetry itself is the consequence of the linear spectrum of excitations and thus 
exists only in the low energy corner. This allows for the use of physically imposed 
deviations from Lorentz invariance, which are governed by the Planck energy scale. 
In a number of papers, therefore, nonlinear dispersion of the quasiparticle spectrum 
has been employed to investigate the quantum physical behaviour related to the 
black hole horizonEm'uEHa. Below, we will find that in the most important cases only 
the first nonlinear correction determines the physics of the horizon, since the largest 
energies involved, although being determined by the Planck energy, are everywhere 
much smaller than the Planck scale. 

Though initially the condensed matter analogue of black hole has been intro- 
duced in a normal liquid, it appeared that the properties of the quantum_yac.- 
uum can be incorporated properly only if one uses quantum supcrfluid liquidaliJ'til. 
The superfluid ground state moving without friction, and an inhomogeneous su- 
pcrfluid velocity v s , provide the quantum vacuum, whose inhomogeneity serves 
as the gravitational field background. The matter propagating in the presence of 
this background is represented by fermionic (in Fermi supcrfluids) or bosonic (in 
Bose superfluids) quasiparticle excitations above the superfluid ground state. They 
form the normal component of the liquid. Two fluid hydrodynamics, introduced 
by Landau and Khalatnikov, incorporates the motion of the superfluid background 
and excitationsEl, and forms the counterpart of the Einstein equations, which in- 
corporate both gravity and matter. A closer equivalence reveals itself when the 
quasiparticles possess "relativistic" properties in their low energy corner, so that 
the superfluid background provides an "acoustic" metrici of Lorentzian signature 
for the low energy excitations. This occurs for example in the case of Bogoliubov 
fermionic quasiparticles in 3 He- A§, and phonons in superfluid 4 He, as well as for 
excitations in dilute Bose condensates. The spectrum deviates from Lorentz in- 
variance for larger energies in a well-defined way, fixed by the exact Bogoliubov 
spectrum in ^e-A, respectively the phonon-maxon-roton spectrum in 4 He. In 3 He- 
A, the analogy is even more attractive, since gravity ariseiL naturally among the 
other low energy collective modes of the sjiperfluid vacuumtil, in the same manner 
as in Sakharov's theory of induced gravityej. 

According to Landau's theory of superfluidity, the motion of the superfluid vac- 
uum is frictionless until a critical velocity is reached, at which quasiparticle creation 
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from the vacuum becomes possible. In supcrfluids with "relativistic"-like spectrum 
of quasiparticles, the Landau critical velocity equals the corresponding "speed of 
light" relevant for quasiparticles, which is the linear slope of the energy spectrum. 
For a generic profile of the superfluid velocity v s -field, "superluminal" motion of the 
superfluid vacuum leads to formation of a horizon for "relativistic" quasiparticles. 
We will call such a horizon "Landau horizon" . If the boundaries of the container are 
situated far enough away from the horizon, the quasiparticle creation predicted by 
Landau for superluminal flow occurs by a quantum process at the Landau horizon 
analogous to the Hawking radiation process at a gravity horizon. 

The value of the Hawking temperature Th = {h/2ir)K s is determined by the 
surface acceleration n s on the black hole horizon. In the black hole analogy of 
fluids!'!, the corresponding relevant quantity is the flow velocity gradient along the 
horizon normal, so that for purely radial flow Th = (h/2Tr)(dv s /dr)\i 10T . In the 
presence of a Landau horizon in superfluids, Hawking radiation leads to quantum 
friction - the superfluid motion is decelerated due to dissipation caused by the 
Hawking radiation processeO'E-i 

The value of the Hawking temperature serves also as a crossover measure to 
distinguish the region of conventional thermodynamics of the quasiparticle gas at 
T 3> Th and the quantum regime T < Th, for which dimensional quantization of 
the quasiparticle energy levels becomes important. The Hawking temperature is 
typically rather small even for condensed matter analogues of the horizon, which 
makes it difficult to observe quantum effects. Therefore, it will prove useful to 
discuss those thermal consequences of the horizon at T > Th which can be simu- 
lated by the thermal states of quasiparticles in the presence of a flowing superfluid 
background. 

In the absence of a Landau horizon the true equilibrium thermal states (i.e. 
those with constant temperature and without entropy production) exist for any ve- 
locity field v s . They correspond to Tolman's law in a gravity field with constant 
Tolman temperature, cf.0, §129. In the presence of the Landau horizon the situa- 
tion changes: True equilibrium states no longer exist. This means that any thermal 
state of superfluid liquid is dissipative in the presence of a horizon. We find such 
dissipative quasi-equilibrium thermal states across the Landau horizon using the 
simplest example of superflow in 1+1 space-time dimensions. They are character- 
ized by constant true temperature T outside the horizon (Tolman's law in gravity), 
and varying T behind it, which satisfies a modified Tolman law. The physics of these 
states is determined by deviations from Lorentz invariance in the ultraviolet region, 
in which quasiparticle excitations can propagate superluminally. The superluminal 
propagation across the horizon allows quasiparticles to exchange information and 
energy, and thus ensures the continuity of the temperature across the horizon. At 
T 3> Th, the dissipation caused by quantum effects related to Hawking radiation 
can be neglected, and the relaxation (and final extinction) of the black hole occurs 
via entropy production behind the horizon, where the "Tolman temperature" is not 
constant. The dissipation is caused by the thermal conductivity and viscosity of 
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the excitation gas (the normal component of the liquid) . 

The high energy corrections to the quasiparticle spectrum are responsible for the 
processes of dissipation in our 1+ld case. These processes are absent in the purely 
relativistic regime of massless excitationsjf the so-called Conformal Killing Vector 
(CKV) conditions are satisfied, see, e.g.^B. These CKV conditions are realized for 
the 1+ld case, where thus dissipation can be caused only by deviations from Lorentz 
invariance in the ultraviolet region. Dissipation is most pronounced in the vicinity of 
the horizon due to the blueshift for quasiparticles there, so that the horizon region is 
responsible for the largest part of the entropy production. This dissipation process 
represents the classical counterpart of the quantum process of entropy production 
by Hawking radiation from the horizon. The same high energy corrections relevant 
for entropy production determine also the entropy of the horizon itself. 

In the second section to follow, we review the properties of the fermionic quasi- 
particle spectrum in 3 He-A and the bosonic one in 4 He. This comprises the influence 
on the energy spectrum of the superfluid background as an effective gravity field, the 
formation of the Landau horizon and the structure of the superfluid vacuum in the 
presence of a horizon, taking into account the nonlinear dispersion at higher energy 
and, finally, the Hawking radiation process in the presence of such a dispersion. 

In the third section, using the "acoustic" Painlcve-Gullstrand metric, we refor- 
mulate the hydrodynamic equations for the normal component in a "rclativistically" 
invariant form, so that they correspond to the thermodynamics of massless 1+1- 
dimensional relativistic particles moving in a gravitational field. 

The thermal states across the horizon are obtained as solutions of the hydrody- 
namic equations on both sides of the horizon in the fourth section, with matching 
between the states occuring due to non-Lorentzian corrections, which are most pro- 
nounced near the Landau horizon. While the equilibrium thermal state outside the 
horizon is determined as usual by a timelike Killing vector, the thermal state be- 
hind the horizon satisfies the general equilibrium condition generated by a spacelike 
conformal Killing vector. 

In section five the entropy related to the thermal states close to the horizon 
is calculated, and dissipation behind and at the horizon is discussed. To compare 
the thermal and quantum processes of the black hole relaxation, we consider the 
quantum correction to the thermal energy in our flow geometry. It has a negative 
sign and thus, as usually assumed, can compensate the diverging term in the stress- 
energy tensor for the Hartle-Hawking state at T = Th- However, we argue that 
this state nevertheless remains dissipative due to the nonlinear dispersion at high 
energy, so that the black hole horizon must eventually shrink at any temperature. 

2. Spectrum and motion of quasiparticles 

2.1. Phonons in a moving superfluid 

A suitably simple framework for our purposes is a superfluid with massless exci- 
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tations, which have a relativistic-like spectrum at low energy. There are two classes 
of such systems: bosonic and fermionic superfluids. A bosonic superfluid is a de- 
generate Bose system which experiences Bosc-Einstein condensation. In the low 
energy limit there are only bosonic excitations in such liquids. In the very low en- 
ergy corner these are phonons, whose energy spectrum is E com — cp, where c is the 
speed of sound, and energy is measured in the reference frame comoving with the 
superfluid condensate. This is realized in superfluid 4 He and in laser-manipulated, 
ultracold Bose-Einstein condensed gases. At higher momentum p the energy spec- 
trum deviates from the relativistic linear behaviour: 

E c ™ = cp + 1 p 3 + ...^cp(l±l (-=2—) +...), (1) 

I Z \-G/Planck/ / 

where 7 can be positive or negative. The parameter 7 is typically determined by 
the analog of the Planck energy -Epianck in superfluids. It will be shown, for positive 
7, that while the first nonlinear correction is important for the horizon problem, the 
higher order corrections are in many cases small and can be neglected. In these cases 
"transPlanckian" physics is not involved, i.e., no "quantum gravity" is necessary. 

The excitation energy in the laboratory frame, in which the superfluid velocity 
v s , the momentum and the energy of quasiparticles are measured, is obtained by a 
Galilean transformation which has the form 

p = p lab = p com , E = E uh {p) =£ com (p)+p-v s , (2) 

or, using ([[]), 

(E-p. v s ) 2 = cV + 2 7 cp 4 + . . . = cV (1 ± (-=2— j I + . . . , :! i 



E P 



lanck 




This equation for the spectrum corresponds to what we would obtain for a plane 
wave mode from a relativistic scalar wave equation in the WKB approximation, 
modified in the comoving frame by higher order derivatives to include nonlinear 
dispersionB. If we neglect nonlinear corrections, we can write the energy dispersion 
(^|) in the form g^p^Pv — 0, where the contravariant components of the metric 
are g 00 = —1, g° l = — u*, g^ = c 2 5 iJ — VgUf. In this Lorentz invariant limit, 
cp -C Ep lanck, the phonons become "relativistic" massless particles of energy E — 
—po, propagating in a (3+l)-dimensional space-time background with metric g^ v 
of Lorentzian signature. This effective metric takes the Painleve-Gullstrand form 

ds 2 = -(1 - {v s /cf)dt 2 - 2c- 2 v s • dvdt + C - 2 (dr) 2 , (4) 

with \J—g = c~ 3 . Observe that we define the metric coefficients such that the 
space-time interval has dimension of time, so that the quasiparticle energy has the 
same dimension as po . 
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2.2. Landau critical velocity, ergosurface and Landau horizon 

If the superfluid velocity exceeds the Landau critical value 

E com (p) 

v L = mm , (5) 

P 

the energy E = E uh (p) of some excitations, as measured in the laboratory frame, 
becomes negative. This allows for excitations to be nucleated from the vacuum. For 
a superfluid velocity field which is stationary in the laboratory frame, the surface 
v s (r) = Vl, which bounds the region where quasiparticles can have negative energy, 
the ergoregion, is called the ergosurface. 

If the dispersion bends upwards, 7 > 0, the Landau critical velocity coincides 
with the "speed of light", vl = c, so that the ergosurface is determined by v s (r) = 
c. In the Lorentz invariant limit, this corresponds to goo( r ) = 0, which is just 
the definition of the ergoregion in gravity. In the case of radial flow (zero axial 
component of v s ) of the superfluid vacuum towards the origin, the ergosurface 
becomes the horizon in the Lorentz invariant limit, and the region inside the horizon 
simulates a black hole for low energy phonons. Strictly speaking this is not a 
true horizon for phonons: Due to the nonlinear dispersion, their group velocity 
v g = dE com /dp = c + 3-fp 2 > c, and thus the high energy quasiparticles are allowed 
to leave the black hole region. It is, hence, a horizon only for quasiparticles living 
exclusively in the very low energy corner, which are not aware of the possibility of 
"superluminal" motion. 

In the case of negative dispersion, 7 < 0, the group velocity v g = dE com /dp < c. 
In such supcrfluids the relativistic ergosurface v s (r) = c does not coincide with 
the true ergosurface, which is determined by v s (r) = vt, < c. In superfluid 4 He, 
where 7 < 0, the Landau velocity is related to the roton part of the spectrum, and 
is about four times less than c.t For inward radial flow, the ergosphere occurs at 
v s (r) — vl < c, while the inner surface v s (r) — c still marks the horizon, in contrast 
to relativistically invariant systems, for which ergosurface and horizon coincide for 
purely radial flow. 

2.3. Painleve-Gullstrand versus Schwarzschild metric in effective gravity 

In our context, it appears worthwhile to point out that in the effective theory 
of gravity, which occurs in condensed matter systems in the form that the primary 
quantity is the contravariant metric tensor g^ describing the energy spectrum, two 
seemingly equivalent representations of the black hole metric, in terms of either the 
Schwarzschild or the Painleve-Gullstrand line elements, are in fact not equivalent 
in terms of the required stability of the underlying superfluid vacuum. Let us 
construct an isotropic Galilean superfluid in 3+1 dimensions, in such a way that 

^The phonon dispersion in superfluid 4 Hc actually first has a slight upturn of phonon dispersion 
(and hence 7 > 0) for small phonon momenta. If we assign 7 < 0, it is thus meant that the 
dispersion after this slight upturn, for intermediate momenta, bends down to the roton minimum. 
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the low energy spectrum of the Bogoliubov quasiparticlcs in the frame comoving 
with the supcrfiuid is given by 

(£ com ) 2 = cV , (6) 

where c is the parameter which plays the part of the speed of light. Then, for 
spherically symmetric velocity profiles v s (r) — =pc^/rg/r, the Painleve-Gullstrand 
line elements are obtained in the laboratory frame as: 

ds 2 = - (l - ^) dt 2 + 2c- l s J^-drdt + c- 2 dr 2 + (r 2 /c 2 ) dQ 2 , (7) 



ds 2 = - (l - — ) dt 2 - 2c- 1 ^drdt + c- 2 dr 2 + (r 2 /c 2 ) dtf 



(8) 



where r§ is the Schwarzschild radius. The Painleve-Gullstrand line elements in 
(@iH) describe black hole and white hole analogues (see e.g. ref.i, a pedagogical 
review of the Painleve-Gullstrand metric is contained i The "surface gravity" 
at the Schwarzschild radius is, therein, k s = dv s /dr\ rs — c/2r$, and the Hawking 
temperature which follows Xh = hK s /2iT. 

An "equivalent" representation of the black or white hole metric is given by the 
Schwarzschild line element, which in terms of the same superfluid velocity reads 

ds 2 = - (1 - v 2 /c 2 ) dP + + (r 2 /c 2 ) dn 2 . (9) 

&■ — 

They are related by the following coordinate transformation: 

t(r, t)=t+( -4t + In ^44) , dt = dt+ -r^dr. (10) 



v s (r) l + v s (r) / 

From the point of view of the effective gravity in superfluids, (0) and (||) describe 
the dynamics of quasiparticles propagating in the superfluid vacuum, moving with 
superfluid velocity =F|u s (r)|, directed to the origin and away from the origin, re- 
spectively. The space-time (t, r) is the absolute space-time of the laboratory frame, 
i.e., as measured by the external experimentalist living in the real world of the 
laboratory. The squared energy spectrum of the quasiparticles, determined by the 
contravariant components of the metric (0,^) , is obtained to be 

(E - v s (r)p r ) 2 = c 2 P 2 , (11) 

or, solving for E, 

E = v s (r)p r ± cp . (12) 

The time t is the time as measured by the "inner" observer at "infinity" (far from 
the hole). The "inner" means that this observer "lives" in the superfluid background 
and makes use of "relativistic" massless quasiparticles (phonons or other excitations) 
to synchronize clocks. The inner observer at some point R ^> 1 sends a quasiparticle 
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pulse at the moment t\ which arrives at point r of t = t\ + f T 1 dr' /\v~\ absolute 
(laboratory) time, where v + and u_ are absolute (laboratory) velocities of radially 
propagating quasiparticles, moving outwards and inwards respectively: 

dr dE , „ , . ,„. 
v± = - = — =±l + v s (r). (13) 
at ap r 

Since from the point of view of the inner observer the speed of light (the propaga- 
tion speed of quasiparticles) is an invariant quantity, and does not depend on the 
direction of propagation, for him or her the moment of arrival of the pulse to r 
is not t but t = (ti + t<i)l% where t<i is the time when the pulse reflected from r 
returns to the observer at R. Since t?, — t\ — f** dr'/\v- 1 + J dr' /\v + | , one obtains 
for the time measured by an inner observer 

t(r,t) = 




(14) 

which is just (|Io|), up to a constant shift. 

At r = rs = 1 one has an event horizon, where the superfluid velocity crosses 
the "speed of light" c = 1. If the fluid moves towards the origin, that is v s (r) < 0, 
this velocity field reproduces the horizon of the black hole (the so-called sonic black 
holeBi): Since the velocity of the fluid behind horizon exceeds the speed c of the 
propagation of quasiparticles with respect to the fluid, the low energy quasiparticles 
are trapped within the horizon. If the fluid moves from the origin outwards, v s (r) > 
0, the velocity field reproduces a white hole horizon. 

In the (necessarily) complete physical space-time of the laboratory, an external 
observer can detect quasiparticles radially propagating into (but not out of) the 
black hole, or out of (but not into) the white hole. The energy spectrum of the 
quasiparticles remains to be well determined both outside and inside the horizon. 
Quasiparticles cross the black hole horizon with the absolute velocity t>_ = — 1— v s = 
—2, i.e. with twice the speed of light: r(t) = 1 — 2(t — to). In the case of a white hole 
horizon one has r(t) = 1 + 2(t — to). In contrast, from the point of view of the inner 
observer the horizon cannot be reached and crossed: the horizon can be approached 
only asymptotically in infinite time: r(t) = 1 + (ro — l)exp(— t). Such a disability 
of the local observer, who lives in the curved world of superfluid vacuum, happens 
because he or she is limited in his or her observations by the "speed of light" , so 
that the coordinate frame he or she uses is seriously limited in the presence of a 
horizon and becomes incomplete. 

The Schwarzschild metric (B) naturally arises for the inner observer, if the 
Painleve-Gullstrand metric is an effective metric for quasiparticles in superfluids, 
but not vice versa. It can in principle arise as an effective metric; however, in 
the presence of a horizon such a metric indicates an instability of the underlying 
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medium. To obtain a line element of the Schwarzschild metric as an effective metric 
for quasiparticles, the quasiparticle energy dispersion in the laboratory frame has 
to be 

£ 2 = c 2 (l^) 2 ^ + c 2 (l-^) P i. (15) 

In the presence of a horizon such a spectrum has sections of the transverse mo- 
mentum p± having E 2 < 0. The consequence of an imaginary frequency of excita- 
tions signals the instability of the superfluid vacuum, if this vacuum exhibits the 
Schwarzschild metric as an effective metric for excitations: Quasiparticle pertur- 
bations may grow exponentially without bound in laboratory (Killing) time, like 
exp[i|9f(2?)|], destroying the superfluid vacuum. Nothing of this kind happens in 
the case of the Painleve-Gullstrand line element, for which the quasiparticle energy 
is real even behind the horizon. The main difference between Painleve-Gullstrand 
and Schwarzschild metrics as effective metrics is thus as follows. The first metric 
leads to the slow process of quasiparticle radiation from the vacuum at the horizon 
(Hawking radiation), whereas the second one indicates a crucial instability of the 
vacuum behind the horizon. 

In general relativity, it is usually assumed that the two metrics can be converted 
into each other by the coordinate transformation in (|Io|). In condensed matter, the 
coordinate transformation leading from one metric to another is not that innocent 
if an event horizon is present. The reason why the physical behaviour implied by 
the choice of metric representation changes drastically, is that the transformation 
between the two line elements, t — ► t + J r dr' v/(c 2 — v 2 ), is singular on the hori- 
zon, and thus can be applied only to a part of space-time. In condensed matter, 
only those effective metrics are physical which are determined everywhere in the 
real physical space-time and come from a physically reasonable energy spectrum 
of quasiparticles. The two representations of the "same" metric cannot, then, be 
strictly equivalent metrics, and we have different classes of equivalence, which can- 
not be transformed to each other by coordinate transformations which are regular 
everywhere. Painleve-Gullstrand metrics for black and white holes are determined 
everywhere, but belong to two different classes. The transition between these two 
metrics occurs via the singular transformation t — > t + 2 J r dr' v/(c? — v 2 ), or via 
the Schwarzschild line element, which is prohibited in condensed matter physics, 
as explained above, since it is pathological in the presence of a horizon. It is not 
determined in the whole space-time, and is singular at the horizon. 

Of importance is the fact that in the effective theory, there is no need for an 
additional extension of the space-time to make it geodesically complete. This is 
because the effective space-time is always incomplete (open), since it exists only in 
the low energy " relativistic" corner, and quasiparticles escape this space-time to 
a nonreiativistic domain when their energy increases beyond the relativistic linear 
approximation regimaiS. 

2.4. Bogoliubov fermions in Fermi superfluids 
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Another class of quasirelativistic systems is provided by Fermi superfluids with 
point nodes in the spectrum, the physically realised example in an earth laboratory 
being the superfluid ^e-A, and possibly chiral p-wave superconductorsE2l. The 
fermionic quasiparticles in 3 He- A are described by a Nambu space Hamiltonian 
which, in the locally comoving frame, has the form 

H A = r 3 e(p) H [np x + T 2 p y ] , e(p) « ^ &v F (p- pf) ■ (16) 

Pf 2m* 

Here, Aa represents the amplitude of the gap, pf and vf — pf/?™* are Fermi 
momentum and Fermi velocity respectively. We assume that the orbital momentum 
vector 1 is oriented along the z axis and that the absolute value of the momentum 
p = |p| is restricted by p — pp <C Pf- 

The Bogoliubov spectrum of fermions in the comoving frame is obtained by 
squaring ( |T(i| ) 

( £co m)2 = £2 = (/ /2m * _4/2m*) 2 + cKpI+pD . (17) 

Here c± = Aa/pf < i>f (the slope of the gap at the node), plays the role of the 
"velocity of light" in the direction perpendicular to 1. 

If the superfluid velocity is restricted to lie in the plane perpendicular to 1 and 
does not depend on z, the momentum projection p z is a conserved quantity, and 
one obtains a 2+ld system with the laboratory frame energies 

ftp \2 (.Pi — Pf) ,22 -2 — 2 2 / -< q \ 

(E-px-Vs) = — 4 ^ 2 +c_ L p_ L , p F =p F -p z . (18) 

In the very low energy regime, i.e., if p is chosen to be close to pf&z, and in 
addition p\ = p^+Py rn^c^, we can write the energy dispersion ( |l8| ) in a Lorentz 
invariant form with (2+l)-dimensional metric (neglecting the quasiparticle motion 
in the z direction). This metric again takes the Painleve-Gullstrand form 

ds 2 = -(1 - (v s /c ± ) 2 )dt 2 - 2cl 2 w s ■ dvdt + c- 2 dr 2 , (19) 

with y/—g = cr 2 . 

The Bogoliubov spectrum, say, for p z = Pf, Pf = reads 

S = v s -px±^/(^-) 2 + cipi. (20) 

Squared, this has the same form as for phonons with positive dispersion in (|J), with 
the "Planck" energy i/pianck = 2m^c 2 l . This shows that the form of spectrum with 
cubic dispersion is generic for any system which has "relativistic" quasiparticles in 
the low energy corner. Further on, we shall consider the positive dispersion, 7 > 0. 

2.5. Physical vacua across the horizon 
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Close to the horizon the horizon surface can be considered to be flat, so that 
the radial flow becomes a one-dimensional flow of superfluid along the normal to 
the horizon (along the x coordinate direction), this flow depending only on x (v s = 
v s (x)x, with v s < 0). Since there is no dependence on y and z, the momentum 
projections p z and p y of quasiparticles are fixed and one obtains (1+1 )-dimensional 
motion of quasiparticles whose "mass" and "speed of light" are determined by p z 
and p y . Figures [l||] display the energy spectrum for the case of zero "mass" , which 
means p y = p z = for the Bose condensate, and p y = 0, p z = ±pp for 3 He-A. We 
measure v s in units of c±; E and p x in units of m*c5_ = Aa(c±/vf) and m*c± = 
Aa/vf, which represent the crossover "Planck energy" and the crossover "Planck 
momentum" to the "nonrelativistic" domain, respectively. Normalised quantities in 
what follows will be designated by a tilde, save for k = p x /m*c± and any velocities, 
which are understood to be scaled by c± (or c in the case of the Bose condensate), 
in the normalised relation: 



E = v s k±^k 2 + ^. (21) 

The position of the horizon is determined by v s {xh) = — 1. We remind the reader 
that, if the nonlinear dispersion of quasiparticles is positive, this is not a true 
horizon, since the group velocity in the comoving frame exceeds the "speed of light" 

1 + — 

,,com| 2 



= -j=±= > 1 , (22) 
'l + f 

and consequently high energy quasiparticles can propagate across the horizon. 

The horizon is understood to be a "Landau horizon" - the position in space 
where the Landau criterium is first violated, i.e., where the superfluid velocity 
reaches the minimum of the comoving frame group velocity: 

v s (x h ) = -mm \v c g om \ = -1 . (23) 

Let us first discuss the vacuum for the Fermi system, which corresponds to fill- 
ing of the negative energy states. The vacuum as viewed in the comoving frame 
corresponds to the filling of states with negative root in ( p0| ) (thick lines in Figj^). 
Outside the horizon, where |v s | < 1 (Fig.|^), the vacuum in the laboratory frame 
coincides with the vacuum in the comoving frame: In both systems the negative en- 
ergy levels are occupied. This, however, does not happen inside the horizon, where 
\v s \ > 1 (FigJ^), and some states, which have negative energy in the comoving frame, 
have positive energy in the laboratory frame. In gravity, the state corresponding 
to the vacuum in the comoving frame is the state empty of particles as seen by 
the "free fall" observer (who moves with superfluid velocity v s ). The vacuum in 
the laboratory frame corresponds to the Boulware vacuum in gravity: It is empty 
for the observer at infinity. In the purely relativistic regime in 1+1 dimension the 
vacuum behind the horizon, determined in the comoving frame, is obtained if the 
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Figure 1: Energy spectrum in comoving frame. The occupied negative energy levels 
designated by thick lines form the comoving frame vacuum. 
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Figure 2: Spectrum in the laboratory frame outside the "horizon" at two points in 
space where v s — —0.5 and v s = —0.8. 



(negative) momentum of upper index, p 1 = (1 — v^)pi — v s p = (1 — vf)p + v s E, is 
used as the Hamiltonian, H = —p 1 , corresponding to the interchange of timelike and 
spacelike coordinates behind a horizon. The states with positive p 1 are occupied in 
this vacuum, and p 1 = — — v^)d x — v s dt] is a conformal Killing vector behind the 
horizon. The same Hamiltonian H = — p 1 determines thermal equilibrium (section 
4.2.). However, any deviation from the linear relativistic regime leads to dissipation 
in this vacuum or the thermal equilibrium state. 

Close to the horizon, the velocity profile can be linearized: v s (x) = — 1 + 
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Figure 3: Spectrum in the laboratory frame inside the "Landau horizon" at two 
points where v s = — 1.2 and v s = —1.5. 



(k s /cx)5x = —1 + x, with Sx = x — x^ <C cj_/k s . Since in the above consider- 
ations we assumed that £ « 1 and x -C 1, only small k are relevant and we can 
approximately write 

fc 3 /8 + kx ~ E . (24) 

The same cubic dispersion has been discussed ini. From equation it follows that 
the characteristic scales for k, E com and x are E 1 / 3 , E 1 / 3 and E 2 / 3 , respectively. 
In dimcnsionful units these are 

E COm = C P ~ 4(a 3 „ck^ 1/3 « Planck , (25) 
\x — Xh\~ — l-g j ~/pianckfy-J f J > ^Planck • (26) 

The characteristic energy, although highly blue-shifted compared with the labora- 
tory energy E, is nevertheless much smaller than the Planck energy scale. Cor- 
respondingly, the characteristic region of the shell (membrane) in the vicinity of 
the horizon, though very narrow, is nevertheless large compared to the Planck 
length scale (the quasiparticle energy in the laboratory frame is determined either 
by the Hawking temperature, or by the temperature of the heat bath, in both cases 
E > Th). Consequently, the higher order corrections to the nonlinear energy spec- 
trum are negligibly small, so that "transPlanckian" physics is not involved save 
for its providing the parameter 7, which can be considered as a phenomenological 
parameter in the effective low energy theory. 

2.6. Wave function 

If the initial state of the superfluid is the vacuum in the comoving frame, then 
in the presence of a horizon this state represents a false vacuum: It possesses large 
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positive energy in the laboratory frame, so that it has to decay. If there is an 
interaction with the walls of the container, the relaxation occurs due to fast processes 
of the dynamics of the superfluid order parameter. However, if the interaction 
with walls is screened, or the wall is very far, the comoving inner observer loses 
the reference frame imposed by the walls. Now the velocity field itself is what 
the inner observer can observe. If the flow is homogeneous, there is no evidence 
for him or her that the superfluid is moving. Thus the only information comes 
from the inhomogeneity of the velocity field, which establishes a special frame, in 
which the velocity profile is time independent. This means that if a horizon is 
present, its relaxation is determined by the gradient of the superfluid velocity. One 
example is quantum Hawking radiation, with the Hawking temperature determined 
by the gradient of the velocity at the horizon. Another example is the classical 
thermal state in the presence of an horizon, whose relaxation is provided by the 
hydrodynamic entropy production induced by the velocity gradient. 

Let us first recall the quantum relaxation. Using the energy E(k) of (^Jj) as a 
Hamiltonian, and canonically quantizing, 



ih , 



ftK s /(m*c^_) , 



(27) 



we write the stationary Schrodinger equation in the momentum representation as 
(x = ihdk): 

iL + nfkdk + ^j i>(k) = Ei>(k) . (28) 

This may be solved by using the ansatz tp(k) — A(k) exp[i$(fc)] to give 

1 



tp(k) = exp 



■ In k + ffT 1 (-E In k + k 3 /24) 



The wave function in i space is, consequently 

1 



ip(x) 



dk exp 



dk exp 



| In k + ifr 1 (kx-E\nk + /fc 3 /24 



1 



In A; + ih 1 [ kx — / x{k)dk 



(29) 

(30) 
(31) 



In the stationary phase (semiclassical WKB) approximation of E ^S> h, the integral 
( |3l| ) will be dominated by three saddle point solutions, which are the three roots 
k(E) of ( pi| ) for given E. Further on, we consider E > 0. For given positive^ 
(outside the horizon), there is only one real solution, k+ s (we use the notation ofu). 
For negative x (inside the horizon) there are three such solutions, fc_ s , k + and 
provided that x is not too small (cf. Fig.[|). The additional ingoing branch in this 
Figure (dotted line), close to the x axis, is obtained if the full initial energy spectrum 
of ( pp[ ) is used. This branch has comparatively high energy for given momentum 
and hence is not relevant for the present discussion. 
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Far enough from the horizon, i.e. at 1 > i > Z? 2 / 3 , we have from ( p4| ) 
k± s ~ E/x, and k± ~ ±2-\/2]5|. The branches k± s — E/x exist in the linear 
dispersion case. At the horizon they approach — oo and +00, respectively. Non- 
linear corrections to the energy spectrum become important when the momentum 
scale E 1 / 3 is reached, and these low momentum branches match with the high en- 
ergy branches k±, which exist due to nonlinearity of the spectrum. The choice of 
the path C in the complex k plane, determining the wave function, was discussed 
in detail iru. 




Figure 4: Quasiparticle trajectories, in the flow v s = —1 + x, for E = 0.1. The 
notations are the same as those iri3: k±(x) are part of the branches with large 
positive and negative momenta and k± s are part of the branches with small positive 
and negative momenta. The arrows indicate the direction of the laboratory frame 
group velocity of the quasiparticle. The horizon is identical with the k axis. Thick 
lines for the quasiparticle trajectory show that the states with momenta k^ s and fc_ 
are occupied in the comoving (free fall) vacuum. The dash-dotted line indicates the 
closest approach distance of the branch of large negative momenta to the horizon, 
at which the negative momentum fc_ part of the wave coming from negative x 
gets Andreev reflected (and the positive momentum k +s part is transmitted). The 
ingoing branch of very small negative momenta is shown by the dotted line. 
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Depletion of the false vacuum occurs by tunneling from the trajectory fc_ • • • fc_ s 
to the outgoing mode trajectory k + ■ ■ ■ fc +s .In the WKB limit of 2-kE/K 3> 1, the 
high energy mode fc_ is Andreev reflectedEll to with a small probability of 
transmission to the outgoing mode. The transmission coefficient is exponentially 
small and corresponds to thermal radiation at the Hawking temperature: 

T ~ exp[-2irE/k] = exp[-E/Tn] , (32) 

where the Hawking temperature is given by the equality Th = h/2ir. The Hawking 
temperature establishes a characteristic scale for the laboratory energy E ~ Th (or 
E ~ Th) and thus the scales for p, E com in (^) and for the thickness of the shell 
xq = \x — Xh\ in J26|). In dimensionless units one has k ~ _E com ~ T^/ 3 , xq ~ T^/ 3 . 

The process of Andreev reflection occurs inside the horizon for our superluminal 
dispersion, whereas it is situated outside for subluminal dispersionB. In both cases 
the Hawking temperature does not depend on the dispersion. The point where the 
wave packet is turning back in co-ordinate space occurs where dx/dk = 0, which is 
given by solving the equation fcp/4 = —Ey/1 + fcf/4. For small fco> there results a 
value fc -(4E) 1 / 3 = -1.59E 1 / 3 , implying x = x(k Q ) = -(i^^+i 2 / 3 /8)E 2 / 3 = 
-0.945i? 2 / 3 . 

2.7. Nonlinear velocity profiles 

An asymptotic form of the velocity profile v s = f(x), where f(x) is some function 
approaching zero (or a constant) at large positive x, and linearity at small Sx around 
the horizon, e.g. f(x) = — 1 + tanh[(55;]), gives us the location of the quasiparticle 
in terms of its wave number as 

x(k) = r 1 (E/k - VTTfcVi) . (33) 

We then obtain a Hamiltonian 

H = E c ™(k) + {kf(x)} , (34) 

the curly brackets specifying that the appropriate operator ordering be taken for 
the expression contained within. 

The influence of a particularly interesting flow profile on quasiparticle motion is 
depicted in Fig.|5|. The profile is uni-directional in the negative x direction, and is 
developing locally superluminal velocities in a depression of the (symmetric) profile. 
It creates a^iharacteristic combination of white hole and black hole (cf. E3o for 
3 He-A and lia for a Bose condensate). Specifically, there are bound quasiparticle 
states inside the horizon, which occur because of the combined effects of Andreev 
reflection, caused by nonlinear dispersion, and the shape of the velocity profile. 
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Figure 5: Large momentum quasiparticle phase space trajectories for the two 
(fixed) energies E = 0.1 and E = 0.5 (dashed), in the velocity profile v s — 
— (3/ cosh 2 [x/a] (/3 > 1), with a = 1, /? = 3. The black and white hole hori- 
zons where v s = — 1 are located at Xh = ±acosh _1 (v / /3) J the upper (lower) sign 
valid for the black hole (BH) and white hole (WH), respectively. The horizon sur- 
face velocity gradient k s = ±(2/a)y/l — The quasiparticle states for negative 
momenta are bound states, where the values of the turning points at large k are in 
leading order k t — ±2/3 for large [3. 

3. Dynamics of the normal component ("matter") 

Now we consider the classical temperature region T 3> Th, where the quantum 
processes of Hawking radiation can be neglected. Still, the behaviour of matter in 
the presence of the horizon is nontrivial and its description requires deviations from 
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Lorentz invariance due to the blueshift in the vicinity of the horizon. 

3.1. "Relativistic" energy-momentum tensor of two fluid hydrodynamics 

The basic conservation equations of s tan dard nonrelativistic Landau two fluid 
hydrodynamics (as expounded in detail ir£3 §8), are those of total momentum and 
energy conservation. Energy and momentum can be exchanged between the two 
subsystems of quasiparticles and superfluid vacuum in a way similar to the exchange 
of energy and momentum between matter and the gravitational field. Moreover, in 
the low energy corner this exchange can be described in general relativistic form. 
The energy-momentum tensor derived from standard conservation equations of two 
fluid hydrodynamics^ can be represented asEU: 

^ = EKf- ^ = - 1+ ?^r^T' (35) 

where / is the quasiparticle distribution function and ^ p = J2 S I d D p/(2irh) D 
designates summation over spin as well as momentum. The group four-velocity is 
defined as 

v* - — v° = l v ■ = 1 dECOm - -(v com Y 

v >° = -( 1 + ?"-^r)' (36) 

Space-Time indices u,v, . . . are throughout assumed to be raised and lowered by 
the Painleve-Gullstrand metric. The group four-velocity is null in the relativistic 
domain of the spectrum only: v^v gfl = if E com = p. Here and in what follows, we 
again set c = cj_ = 1 (scale velocities with c), except where the "speed of light" is 
explicitly indicated for clarity. 

The energy-momentum tensor thus defined satisfies the covariant conservation 
law equation 

As pointed out by Landau and Lifshitz@, this equation does not generally express 
any conservation law whatever. The same is true for quasiparticles in superfluids. 
This equation does not mean that energy and momentum of quasiparticles are 
necessarily conserved: In the presence of an inhomogeneous condensate the second 
term describes the energy and momentum exchange between quasiparticles (matter) 
and superfluid vacuum (gravity field) 

{V—gT* 1 ,;)^ = P t d v v l s + V fd„E com , (38) 



p 

where P = X) P P/( r >P) ^ s ^ nc quasiparticle momentum density. The last term 
appears only if the "speed of light" depends on space-time. Further on we neglect 
this term. 
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In a onc-dimcnsional situation, and assuming that there is no temporal depen- 
dence of superfluid velocity, respectively if it is permitted that it is neglected, we 
obtain an energy conservation law for quasiparticles in the form 

|(V^T° ) + ^(^r x o)-0, (39) 

which tells us that in any time independent situation, the energy flux ^J~ r g~T x a is 
necessarily constant in space. The momentum conservation in x direction reads 

^- t (V=gT° x ) + ^-{J—gT x x ) = -P x 8 x v s = -±T° x V=g. (40) 

It contains the velocity gradient as a source term, which corresponds to the "gravity" 
k = d x (vl/2), where plays the role of the gravitational potential. 

The relevant^ components of the energy-momentum tensor are (po = —E, p° = 

E com ): 

\/—gT Q i = ~^^.{Pi = Pi Momentum density in either frame, 
p 

gT° = - y^.fE (Negative) Laboratory frame energy density, 
p 

gT k i = fPiVg Laboratory frame momentum flux, 
p 

—gT l = = - ^ fEv l q (Negative) Laboratory frame energy flux, 
p 

— T oo = ^fp = '^2fE com Comoving frame energy density. (41) 
p p 

We will further on keep the square root of the negative metric determinant 
for structural clarity, though it is a space-time constant in our considerations, as 
we assume the "speed of light" to be constant below. The relation between the 
laboratory frame energy density and that in the comoving frame expresses the 
underlying Galilean invariance: 

-V=gT° = f E = ^~9T m + vlPi . (42) 
p 

The additional hydrodynamic equation to be satisfied is that of the conservation 
of the number of particles comprising the underlying liquid, say, the number of 3 He 
atoms. If temporal variations of the density are negligible, this amounts to requiring 
that the total momentum density 

j = pv s + P (43) 

i "Relevant" are termed those components of the energy-momentum tensor that do not depend on 
the low energy properties of the quasiparticles, i.e. on the form and existence of a Lorentz metric. 
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has vanishing divergence, V • j =0. That we can neglect (temporal and spatial) 
variations of the density will be true in the case of 3 He-A, where the "speed of 
light" cj_ is much smaller than the sound velocity. This is in marked contrast to 
truly "acoustic" black holes, as they will exhibit large variations of the density for 
velocities of the fluid above the speed of sound. 

3.2. Components of energy-momentum in local thermodynamic equilib- 
rium 

Local thermodynamic equilibrium is characterized by the local temperature T 
and local velocity v„ of the quasiparticle system. The equilibrium distribution 
function isEH: 

1 1 

fr = = (44) 

J l + cxp[T- 1 (£;(p)-p-v„] 1 + cxp[T- 1 (£; com (p) - Pi w l {x))} v ; 

where w = v n — v s is the so-called counterflow velocity. 

In local thermodynamic equilibrium the components of energy-momentum for 
the quasiparticle system (matter) are determined by the generic thermodynamic 
potential (the pressure), which for fermions has the form 

n = - T ^E/^Ml-/). (45) 
The momentum density is expressed in terms of Q as 

= v^n = i£^ p/ft= «U {Pnhw * . (46) 



Here p n is the so-called density of the normal component (in general a tensorial 
quantity), so that the total current of the liquid: j = pv s +p n (v n —v s ) = 
where p s = p — p u is the density of the superfluid component. 

Considering below for simplicity the (spatially) one-dimensional case (D = 1), 
one has for the purely "relativistic" contribution such that E com = \p x \ (and thus 



"-Ess*. r " = 7rb»' m 

2fi 



Px = V^gT° x = p n w , p n = 2 . (48) 

The other components are in the "relativistic" case then given as 

= (e + tt)rt l u u + n g " v , e = -fi + T— = DVL, T% = . (49) 

where u a and u a — g a pu^, which satisfy the normalisation equation u a u a = — 1, 
are 



1 i V) w i l + w-v w 

Vl — w 2 v 1 — w 2 yl — w 2 v 1 — w 2 



(50) 
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The distribution of quasiparticles in local equilibrium in ( f44| ) is expressed via the 
temperature four- vector [3^: 

For the relativistic system, true equilibrium with vanishing entropy production is 
established if (3^ is a timelike Killing vector: 

Pvn» + = . or P a d a g^ + {g„ a d v + g V ad„)(3 a = . (52) 

For a time independent, space dependent situation the condition = po-,o = (3 l digoo 
gives /3* = 0, while the other conditions are satisfied when (3° = constant. Hence 
true equilibrium requires that v n = in the frame of the velocity texture, and 
T = constant. These are just the equilibrium conditions in superffuids. They 
correspond to Tolman's law in general relativity, with Tolman temperature being 
equal to the real temperature of the liquid: 

T T 

T ofr = = -== . (53) 



4. Thermal states in the presence of a horizon 



4.1. One- dimensional hydrodynamic equations 

To describe the thermal states in the presence of the horizon in a most sim- 
ple way, we consider a purely one-dimensional situation. This is, however, not 
completely unrealistic an assumption. Of course, we have to relax the condition 
j =const., which does not allow for a velocity v s gradient in the one-dimensional 
case. This, however, can be achieved by the inhomogeneous profile of a tube. Then 
there is some flow in the transverse directions, such that V ■ j = is satisfied. 
An alternative possibility to construct a one-dimensional horizon is to make cj_ (x) 
space dependentliJ. The mass conservation law j = const, also implies that the 
nonzero quasiparticle momentum, which arises, say, due to Hawking radiation, has 
to be compensated by a corresponding change of the superflow velocity. This is 
the mechanism of back-reaction of the superfluid vacuum to quasiparticle dynam- 
ics. Back-reaction is, however, small and can be neglected, namely if the effective 
temperature of quasiparticles is small compared to the superfluid transition temper- 
ature T c . This will always be true in our approximation, for which T < m*c^ <C T c . 
Then, P < pv s . 

We will investigate below a case for which the quantum effects related to gravity, 
including the Hawking radiation process, can be neglected. This is permitted if all 
the relevant energies are much higher than the Hawking temperature: T ^> Tjj, and 
H/t Tjj. In the latter relation, / = cr is the mean free path, and the relation 
means that I is small compared with the characteristic length, within which the 
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velocity (the "gravitational potential") changes: l(dv s /dx) <C 1. This is just the 
condition for the applicability of the two fluid hydrodynamic equations, where the 
variables are the superfluid velocity v s which, when squared, plays the part of 
the gravitational potential, as well as temperature T(x) and velocity v n (x) of the 
normal component, which characterize the local equilibrium states of quasiparticles 
(matter). The dissipative terms in the two fluid equations can then be neglected 
in a zeroth order approximation, since they are small compared to the reversibly 
hydrodynamic terms by the above parameter l(dv s /dx) -C 1. 

If the superfluid velocity (gravity) field is fixed, the other hydrodynamic vari- 
ables, temperature T(x) and velocity v n (x) of "matter", are determined by the con- 



servation of energy and momentum. From equation (39) it follows that the energy 
flux Q carried by the quasiparticles is constant. In the relativistic approximation 
one then has 

t v n (l + wv s ) _ _ ^ _ (2s + l)^ 



Q = -V^9T X = 2ft^- = const , ft = v ' ' T ff . (54) 

I — ur 12 

From ( [Io|) there results the first order differential equation 

/ f\ fn \ /in 

-d x [2Q + ft = 2ft- ^d x v s . (55) 

V 1 — «r I I — w z 



4.2. Local equilibrium states with zero energy flux 

In the case that the energy flux is zero, equation (^34|) yields two possible states. 

4.2.1. True equilibrium outside horizon 

Equation (|54j) is satisfied by the trivial solution v n = 0. Then from ( p5| ) it 
follows that T = constant. This corresponds to a true equilibrium state outside 
the horizon. The effective temperature satisfies Tolman's law in (f^): T c s(x) = 
T/y— goo(x) = T / y/l — v%(x). This equilibrium state cannot be continued across 
the horizon: The effective temperature T c g diverges when the horizon is approached 
and becomes imaginary inside the horizon, where \w\ — \v s \ > 1. 

4.2.2. Thermal state behind horizon 

Equation ( |54| ) is satisfied by I + wv s = (or u Q = 0). Since w 2 < 1, this 
solution can be valid only inside the horizon, where v 2 > 1. From j5q) it follows 
that ft = const/ (v 2 — I), and thus the temperature behaves as T 2 oc ft(l — w 2 ) = 
const/fg , or T — Th or) i e ft/|« s |, where Th or ,ieft is the temperature at the horizon when 
approached from inside. Thus inside the horizon one has a quasi-equilibrium state 
with inhomogeneous temperature. The effective temperature behind the horizon 
follows a modified Tolman law: 
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The above equation is valid only in 1+1 dimensions. For larger spatial dimension 
D > 1 one has Q oc {v 2 - \)-(d+i)/2D and T ^ K ^2 _ 1 )-i/2£> ) so that T K 

The thermal state behind the horizon is no longer determined by the Hamiltonian 
H(x > 0) = — po = E, but can be described as a local equilibrium state if one uses 
the kinetic momentum as the Hamiltonian: 

/(p, x < 0) = — t n w^ r . ( 57 ) 

1 + exp(H(x < 0)/r hor joft) 

H{x<0) = ~ P \ p 1 = {1 - v 2 s )p + v s E = p + v s \p\ = -p°x . (58) 

Here x = dE/dp, so that p 1 is a "mass" multiplied by a velocity and thus represents 
a kinetic momentum. The reason why the kinetic momentum plays the part of the 
Hamiltonian is that 

_p!_=/3M( x< o)p„ ) p=(-^,±ZlL) . (59) 

J hor,left \ J hor,left J hor,lcft/ 

Here (3^{x < 0) is a two-temperature, which satisfies the CKV condition + 
Pu;u ^ 9nv It can be verihed by direct calculation that 

+ Pu;n = 9^9 X {1 - v 2 ) . (60) 

For massless relativistic particles the CKV condition is sufficient for thermal equi- 
librium states without dissipation to exist. Thus, despite the inhomogeneous tem- 
perature and nonzero normal velocity v n , the state behind the horizon becomes a 
true equilibrium state in the low energy relativistic limit. 

Existence of nondissipative thermal states behind the horizon is a consequence 
of the 1+1-dimensional situation, in which the metric is conformally flat. The co- 
ordinate transformation i = t — J x dx'v s /(l — v 2 ) , x = x transforms the vector 
/3 M — > /3 M = (0, 1 — v 2 ), so that p 1 = — i(l — v 2 )dx- Another coordinate trans- 
formation dx/(l — v 2 ) = dx transforms p 1 to —id$. On the other hand, these two 
transformations give the conformally flat metric ds 2 — (l — v 2 )(—dt 2 + dx 2 ) with the 
conformal factor 1 — v 2 . This implies that p 1 is a conformal Killing vector satisfying 
the CKV condition in (pOj). 

4.2.3. Matching thermal states across horizon 

Since the superluminal dispersion provides an energy exchange between the mat- 
ter (quasiparticles) inside and outside the horizon, the temperature must be con- 
tinuous across the horizon. Thus Thor,icft = ?hor, right- It is important to note that 
even in the comoving frame, the blueshifted energies are finite and much smaller 
than the "Planck" scale, on which the hydrodynamic formalism breaks down. Con- 
sequently, the temperature and the normal fluid velocity can be well determined 
not only on both sides of the horizon but even at the horizon itself, where all the 
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characteristic scales are determined by the first nonlinear correction to the quasi- 
particle spectrum. In our 1+ld case thermal states on both sides of the horizon do, 
however, not depend on the details of this thin layer. The only effective output of 
the membrane is that it establishes the boundary condition Th or ,ioft = ?hor, right- 

Thus, assuming that at +00 the superfluid velocity vanishes, one has the follow- 
ing temperature profile on both sides of the horizon: T(x > 0) = X^, T(x < 0) = 
Too/\v s (x)\. The distribution functions on both sides of the horizon have the form 

x>0)= i + cxpHVToo) ' Po = W + VsP ' (61) 

/(p,.<0)= 1 + exp( 1 _ pl/Too) , pi=v a \ P \+P. (62) 

There is a symmetry between the states inside and outside the horizon. While 
outside the horizon one has Tolman's law, T(x > 0) = T fF J I — v 2 , the temperature 
behind the horizon is obtained from Tolman's law by the substitution v s — > l/v s , 
i.e., T(x < 0) = T e g^/1 — The same procedure applies to w: One has w = —v a 

outside the horizon and w = —l/v a inside. The temperature two- vector for the state 
behind horizon 

Fix < 0) = = ^(-v s , 1 - v 2 a ) , ^{x < 0) = (/3 , 0i) - ^-(0, 1) , 

(63) 

can be compared with the temperature two-vector for the true equilibrium outside 
the horizon 

/n*>0) = i(l,0), ^(x>0) = ^(v 2 s -l,-v s ) . (64) 

The relation between these two-temperatures shows that the states across the hori- 
zon are dual to each other: 

^(x > 0) = e^p^x < 0) , (65) 

i.e., the space and time components of the two-temperature transform to each other 
across the horizon (e 11 = e 22 = , e 12 = — e 21 = 1). It also appears that the thermal 
equilibrium behind the horizon is the same as the thermal equilibrium outside the 
horizon, if the superfluid velocity and the speed of light are interchanged. Behind 
the horizon the speed of light plays the part of the superfluid velocity and vice 
versa. This corresponds to the case of a soliton moving in superfluid 3 He-A, where 
the speed of light changes in space, while the superfluid velocity remains constants. 

4.3. Local equilibrium states with nonzero energy flux 

If the energy flux is nonzero, we express f2(u s ,iu) = Q(l — w 2 ) /2v n (l + wv s ) 



using (p4[), and substitute it into fl55|). One obtains an equation for w(v s 
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which has the solution 

(1 - v s ) - C(l + v s ) 
where C is a constant. Using this equation we have the sequence of relations: 
v n = w + v s = (l-C)- 
l + wv s = (1 + C)- 



' (1 - v s ) + C{1 + v s ) 
1-v* 



{l~v s ) + C(l + v s ) 
= (l-tfi^. (68) 

Since \w\ < 1, one has C+ > outside the horizon, where \v s \ < 1, and C_ < 
inside the horizon, where |t> s | > 1. The energy flux is 

Since the energy flux can be only inward, Q < 0, (in the absence of Hawking 
radiation), one has C 2 > 1. From (69) it follows that the effective temperature 
exhibits Tolman's law'*, determined by the temperature at +oo (if u s (+oo) = is 
assumed) : 

T cS = 2| , Tl = \Q\ - 24 ' C| . (70) 

<y\l-v%\ (2s + l)7rC 2 -l 

It follows from the above relation that a modified form of Tolman's law, telling 
us that the "temperature" detected by a local observer, T cff = const./ y/\ — gool = 
const. — vf\, is fulfilled even in the absence of complete equilibrium. "Temper- 
ature" is set in quotation marks, because T eS represents only an effective "relativis- 
tic" temperature of the quasiparticles appearing in the distribution function, while 
the true temperature of the superfluid still remains T. The effective temperature 
diverges at the horizon, while the true temperature remains finite, but changes in 
space: 



T= VI - w A T cS ^Too = Too -r: , -7, 71 

V|l-«f] \l-v B + C(l+v s )\ 

and approaches at the horizon the value 

Thor^TooVlq. (72) 



Since the true temperature must be continuous across the horizon, the condition 
C_ = — C+ is enforced. 



§The relativistic Tolman law of a Lorentzian space-time obviously involves corrections from the 
nonrelativistic tail. These corrections are, however, small far enough away from the horizon. 
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Hence the continuity of T and the conservation of the energy flux completely 
determine the thermal state inside the horizon in terms of the temperature at posi- 
tive infinity, Too, an d the energy flux Q. In terms of Too and t; n (+oo), the relations 
v n (+oo) = (1 - C+)/(l + C+) or C+ = (1 - w n (+oo))/(l + u n (+oo)) hold. In the 
limit of vanishing flux, \C\ — > 1, one obtains the results of the previous subsection: 

outside horizon : T = Too , v a = , w(x) — —v s (x) , (73) 

Too v 2 (x) — 1 1 

inside horizon : T(x) = °° , v n {x) = — — -—— , w(x) = — . (74) 

\v s (x)\ v s (x) v s {x) 



5. Properties of the horizon 

5.1. Kinks at the horizon 

In the presence of superluminal dispersion, all physical quantities are contin- 
uous at the horizon. On the other hand, the thickness of the membrane, deter- 
mined by the nonlinear disperion in (|2^) with E = T^, is very small and can 
be neglected. In this limit the thermal quantities experience kinks. For exam- 
ple, VT|+o = VT|_ (C+ - 1)/(C_ - 1); Vw| ±0 = -C±Vu s (x = 0). In the case 
of zero flux, equations ( |73"| - |7"4| ) , the jump in the derivative of the temperature is 
VT| + o — VT|_o = ^ttTqqTh, and does not depend on details of the dispersion. This 
means that in the limit of a purely relativistic system, the presence of a nonzero 
temperature at infinity implies a singularity at the horizon. This singularity at 
the horizon cannot be removed, since in the presence of "matter" with nonzero 
temperature the system is not invariant under coordinate transformations, and this 
produces a kink in temperature at the horizon. 

5.2. Entropy related to horizon 

The thermal states also have an entropy profile, for which we can distinguish the 
part coming from the bulk and that from the horizon. Typically the entropy related 
to the horizon is dominating, since it involves the "Planck" scale parameter. We 
again discuss the 1+ld case, where the horizon contribution is only logarithmically 
larger than the contribution of the bulk. 

The true entropy of superfluid, which is carried by its normal subsystem at 
^ » T H , is 

S = Jd D rS, S = ^. (75) 

The "relativistic" entropy, which is measured by a local observer living in the quasi- 
particle world is the effective entropy 



(76) 
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In the presence of a horizon, the total true entropy can be divided into 3 contribu- 
tions: 

S = S cx t + <Si n t + 5hor ■ (77) 

The exterior entropy, which comes from the bulk liquid, is proportional to the size 
L oxt of the external region: iS ex t oc T^L^. A similar estimate holds for the entropy 
of the interior region, S- ln t oc T^L^t- The entropy related to the horizon stems from 
the logarithmically divergent contribution at the horizon: 

5 hor = ^t 1 ^ ( r° + I 1 '"*) dsTT^M 1 + (78) 



6 \J-i/ Ks Jx j ii-f s 2 r 

where we have used relations (|73|) and (|74|) (valid for the limes of vanishing energy 
flux). According to (|24|), the ultraviolet cutoff xq is provided by the thickness of the 
shell (membrane) determined by the nonlinear dispersion, which gives (see equation 
@ with E = T ao ): 

x a n s =i ~ f^ 3 = ' ■ (79) 

\ ^Planck / 

5W=£^§Sm(^g4) . (80) 



As a result one has 



This relation also means that the density of quasiparticle states diverges logarith- 
mically at the horizon: 

= < 2 » + " / ^ - - W - • 

T H « E « £ Planck . (81) 

5.3. Dissipation at the horizon 

Behind the horizon the conditions for the true equilibrium in ([52] ) are not satis- 
fied, it M „ = I3 U - V + (3 u;IJ , 7^ 0, cf. equation (|dC|): 

u 00 = ^(l-v*)d x (v 2 s ), "oi = ^d x (v 2 s ) , uu = ~d x tf) . (82) 

But quasi-equilibrium CKV conditions are met, with the conformal function being 
equal to the local gravitational field, normalized to the temperature at infinity: 

u^ = 2g^, $ = -^-d* (^) • (83) 

Since the true temperature and the normal component velocity, T and v n , are space 
dependent, there is entropy production and dissipation stemming from thermal 
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conductivity, giving a term ^(VT) 2 , and second viscosity, contributing a term 
C2(divv n ) 2 . In the truly relativistic regime with massless fermions dissipation is, 
however, absent since the CKV conditions are satisfied. This fact isalso known 
for the bosonic case, superfluid 4 He, where, according to KhalatnikovE 2 ], these two 
transport parameters are zero in a pure phonon gas. They become nonzero for 
any deviation from the linear phonon spectrum. Thus, if the Hawking process is 
neglected, in the purely relativistic regime, there will be global equilibrium states 
on both sides of the horizon, for any given temperature at infinity. Dissipation 
occurs only due to deviations from the relativistic spectrum. These deviations are 
most pronounced and effective within the thin shell of the horizon vicinity, so that 
the main source of dissipation in a 1+ld superfluid, in the presence of a Landau 
horizon, will be concentrated at the membrane. The dissipation is determined by 
the gradient of the superfluid velocity at the horizon. 

This, then, is the thermal counterpart of Hawking radiation. The latter becomes 
important only when the quantum limit is approached, i.e., when approaches 
the Hawking temperature Xh- The hydrodynamic dissipation at T x 3> Xh does the 
same job as Hawking radiation at X = 0, leading to the final extinction of the black 
hole. In our case, this means deceleration of the superfluid below |k s | = 1 and thus 
the consequent merging of black and white hole horizons (cf. Fig.^j) , and, finally, 
their mutual annihilation. 



5.4. Quantum corrections 

When the external temperature is low, and comparable to the Hawking tem- 
perature Xh, the quantum correction to the effective action becomes important. 
The quantum part of the action, representing outside the horizon theenergy of the 
Boulware vacuum, is given by (see for example equation (30) in RefE3): 

Quantum = W^^^ft^ ' ^ 

where R is the Ricci scalar. If we are interested in the time independent action, 
then the metric g^ v and the curvature R are time independent as well. Let us first 
consider the space outside the horizon, where we have 

d^"d v =d x {l-vl)d x . (85) 

After an integration by parts, we obtain 

lWum(x > 0) = - _L_ ^ ^ R ) (fl-lfl) , ( 86 ) 



1927T I - w s 2 

The curvature R calculated by using the "acoustic" metric outside the horizon yields 



R = d 2 x {v 2 s ), 



(87) 
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2S +1 ,„ , 2 2 1 



and thus the quantum correction becomes 

2s 

^quantum(x > 0) = (d x {vj ) ) - 

This coincides with the result of Ref@, obtained outside the horizon for the Boul- 
ware vacuum. Like in the case of the thermal energy, we may extend this result 
to the region behind the horizon. Behind the horizon, one has to consider the vac- 
uum state as being the limit of vanishing normal component contribution (the heat 
bath) , which keeps the value of the normal component velocity v n fixed. This means 
that we have to consider the vacuum state with respect to the moving heat bath. 
As a result, the coordinate dependence of the quantum action behind the horizon 
is determined by the counterflow velocity w, instead of the superfluid velocity v s . 
Since w = — l/v s (in the limit of vanishing energy flux), the result is given by (|88|), 
using the substitution — > 

O quantum (x < 0) = . (89 ) 

Close to the horizon both regions, outside and inside the horizon, therefore con- 
tribute in a similar way: 

"quantum ~ ^ h_ ^21 ' ^ > 

The total energy density close to the horizon, quantum + thermal, is given by 

O(total) = ^(quantum) + fi(T) m i — 2_J±^2 Ji (91) 

12 |l-f s 2 | 

In agreement with more general theory (see, e.<7-,0), the diverging terms in the 
thermal and quantum energy compensate each other if the external temperature at 
spatial infinity equals the Hawking temperature, = Th, that is, in the Hartlc- 
Hawking state (the stress-energy tensor in the Hartle-Hawking state, defined by 
the vacuum of the freely falling observer, is regular at the horizon). However, 
despite such a compensation, the Hartle-Hawking state remains dissipative: The 
true temperature T is inhomogeneous behind the horizon and the velocity v n of the 
normal component is nonzero. This leads to dissipation behind the horizon, due to 
the high frequency nonlinear dispersion. 

The quantum correction does not change the entropy of the horizon in the ther- 
mal state of (|S0|). But in the particular state of = Th, the Hartle-Hawking 
state, the quantum correction ancLthfijant ropy of the horizon are related by a gen- 
eral theorem. According to Refs.EaOO, this entropy is expressed in terms of the 
Lagrangian as 

f dL 
Shor = -87T / t u n v txn p — da , (92) 
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In the 1 + 1-dimcnsional case, for which R = 2i?oioi/ff = — 2i?oioi> it follows that 



Shor = — 87T- 



9-Rqioi 

Using for the Lagrangian, one obtains 



= 16tt 



hor 



dR 



hor 



She 



(2a + 1) 
12 
(2s + 1) 



1 



|1-^ S 2 I x 
(ln|l-, s 2 |) hor = 



hor 



., (2s + 1) -Epianck 

~ 11 - m ■ 



(93) 



(94) 



12 v 1 sl/hor 18 T H 

where we took into account the cutoff induced by the high-energy dispersion, |1 — 
Ug| ~ (JH/i^pianck) 2 ^ 3 - The expression (94) coincides with the estimation of the 
entropy of the horizon in equation (pC|), if T^ = Th, i.e., for the Hartle-Hawking 
state. 



6. Concluding Remarks 

Horizons and ergpregions can be simulated in a multitude of different systems 
(cf. the overview inES). The most promising systems, displaying quantum behaviour 
related to a nontrivial effective space-time, are superfluid 4 He, superfluid "Tie, and 
dilute Bose-Einstein condcnsatesE^tJ, which can be conveniently manipulated by 
lasersB 

The Hawking-analogous process itself is exceedingly difficult to measure in exper- 
imental practice. The Hawking temperature given by 2ttTh = hn s is, in experimen- 
tally relevant units, Th = 1.22 nKXK s [10 sec ] . In 3 Hc-A, where the velocity gradi- 
ent attainable, in principle, is limited by k s < c_l/£a> this implies that Th < 0.1 /iK. 
In superfluid 4 He, the corresponding limitation, k s < c/£, for Th is much less severe, 
and comparatively large Th would be available (Th ~ O(10 mK) ~ O(10~ 2 T c )), for 
k s ~ O(10 _2 c/^). The Hawking process would, however, probably be masked by 
non-Hawkingian density perturbation excitations, generated by trying to force this 
dense fluid to move with large gradient of velocities near the speed of sound. Mat- 
ters will possibly improve for a dilute Bose-Einstein condensate, e.g. a sodium 
vapourEI, T H - 0(0.1 nK) - O(10~ 4 T c ), for k s - O(10- 2 c/£), because there to 
sustain this large flow gradient without generating other perturbations not related 
to the Hawking process might be simpler. 

It appears, however, that the classical dissipative signature of a horizon for 
temperatures well above Th, which we have been discussing in section 5.3, will 
more straightforwardly be detected with current technological means. The results 
of the present work can be extended to the more realistic 2+1 and 3+1 cases. In all 
cases the thin shell, which includes both sides of the horizon, is of great importance. 
In the shell the large blueshift is compensated by the first nonlinear correction to 
the linear "relativistic" energy spectrum. The thickness of the shell is much smaller 
than the scale of the velocity profile (the radius of the horizon in the black hole 
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analogy), but also much larger than the microscopic length scale (the Planck length 
in the analogy) . This allows for a determination of the temperature of the thermal 
state everywhere across the horizon. This shell is responsible for the most important 
properties of the thermal states in the presence of the horizon, including the entropy 
and entropy production, both concentrated in or close to the shell. 

In principle even the 1+1-dimcnsional case discussed here can be realized. The 
condition for the applicability of our 1+1-dimcnsional results is that the transverse 
energy level separation for, say, a torus geometry of the BEC be larger than the 
temperature, such that these levels are not (appreciably) occupied. Such a torus 
trap may be realized by using a highly anisotropic trap, whose trapping frequency is 
"soft" in only one, ring-shaped, direction. The condition on the transverse energy 
level separation reads T <C he/ '2d, where d is the diameter of the torus of cross 
section a = 7rd 2 /4. For that same sodium BEC we have been using for an esti- 
mation above, for c/2d ~ 10 3 sec _1 (c ~ 0(1 — 10 mm/sec), d ~ 0(/im)), it turns 
out that T < 10 nK ~ 0(0.01 • • - 0.1T C ) is to be fulfilled, which appears feasible 
(recently improved calorimetry has already made it possibles to measure temper- 
ature changes of less than 10 nK). A closer investigation of this classical dissipativc 
process in these systems should, under general conditions, then take into account 
the effect of back-reaction on the phenomenon of horizon entropy production. 
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